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ABSTRACT 

The detection of astrophysical very high energy (VHE) neutrinos in the range of 
TeV-PeV energies by the IceCube observatory has opened a new season in high energy 
astrophysics. Energies ~PeV imply that the neutrinos are originated from sources 
where cosmic rays (CRs) can be accelerated up to ~ lO^^eV. Recently, we have shown 
that the observed TeV gamma-rays from radio-galaxies may have a hadronic origin 
in their nuclear region and in such a case this could lead to neutrino production. In 
this paper we show that relativistic protons accelerated by magnetic reconnection in 
the core region of these sources may produce HE neutrinos via the decay of charged 
pions produced by photo-meson process. We have also calculated the diffuse flux of 
HE neutrinos and found that it can be associated to the IceCube data. 

Key words: Low luminosity AGNs- very high energy neutrinos- cosmic ray acceler¬ 
ation: magnetic reconnection- radiation mechanisms: non-thermal. 


1 INTRODUCTION 

Neutrino observations can provide unique information to un¬ 
derstand their origin and can even lead to the discovery of 
new classes of astrophysical sources. The inherent isotropic 
nature of the detected neutrino flux by IceCube is compati¬ 
ble with an extragalactic origin and is supported by diffuse 


high energy y-ray data (Ahlers & Murase 20141. The ob¬ 


served neutrinos with energies ~ PeV suggest that they are 
originated from a source where cosmic rays (CRs) can be 
accelerated up to ~ lO^^eV. 

A potential mechanism to produce VHE neutrinos in 
the TeV-PeV range is through the decay of charged pions 
created in proton-proton (pp) or proton-photon (py) col¬ 
lisions in a variety of astrophysical sources which, in the 
framework of the IceCube observations, may include active 


galactic nuclei (AGNs) (|Kasanas 

& Ellison||1986 Stecker et 

al.|1991 Atoyan & Dermer|2001 

Neronov & Semikoz|2002 1 

and gamma-ray bursts (GRBs) ( 

Waxman & Bahcall|1997|). 


Hadronic mechanisms producing VHE neutrinos via the 
acceleration of cosmic rays (CRs) in AGNs have been sug- 


gested for more than three decades (|Eichler|1979||Protheroe 

& Kasanas |1983 Mannheim| 1995| Hazlen & Zas 

1997 

Mucke & Protheroel 20011 Kalashev et al. 20141 Marinellil 

& Fraija|2014b 

Atoyan & Dermer|2003 

Becker|2008 

. Gur- 
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rently, the detection of gamma-ray emission at TeV energies 
in AGNs, not only in high luminous blazars, but also in less 
luminous radio-galaxies, has strengthened the notion that 
they may be excellent cosmic ray accelerators and therefore, 
important potential neutrino emission candidates. 

Several recent models have tried to describe the de¬ 
tected TeV neutrino emission as due to AGNs. For instance, 


Marinelli & Fraija (2014b I employed two different hadronic 


scenarios involving the interaction of accelerated protons at 
the AGN jet either with photons produced via synchrotron 
self-Compton (SSC) or with thermal particles in the giant 
lobes. They then derived the expected neutrino flux for low 
luminous AGNs (LLAGNs) Q or more specifically, for radio 
galaxies for which they examined the origin of the observed 
TeV gamma-ray spectra as due to hadronic processes. 


Earlier work by Gupta (20081 had already introduced 


hadronic scenarios to explain the TeV emission in LLAGNs 
(e.g., Cen A). Also, [Waija| ( |2014a|b[ ) suggested neutral pion 
decays from pp and py interactions in these sources as prob¬ 
able candidates to explain the high energy neutrinos. In an¬ 


other model, Kalashev et al. (20141 attempt to reproduce the 


^ By LLAGNs we mean non-blazar sources with Lrr^ < 
10‘^°ergs“l (see Ho et al.|l997 Nagar et al.|2005l, where is 


the Ha luminosity. These typically consist of liners and seyfert 
galaxies which are also FR I or FR II radio sources. For more 


details see Kadowaki et al. (2015i. 
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IceCube data using the P 7 mechanism considering the radi¬ 
ation field produced by the accretion disk around the AGN 
central black hole (assuming a standard Shakura-Sunyaev 
accretion disk model). Alternatively, Kimura et al. (20141 
calculated the neutrino spectra using the radiatively ineffi¬ 
cient accretion flows (RIAF) model in the nuclei of LLAGNs 
considering both pp and py mechanisms and stochastic pro¬ 
ton acceleration in the RIAF turbulence. 

The possibility of producing VHE neutrino emission has 
been also extensively explored in blazars - AGNs for which 


the relativistic jet points to the line of sight (e.g., Atoyan 


& Dermer 

|2003| |Becker|2008| IMurase et al.|2014||Dermer et 

al. 

2014 

I. 

Dermer et al. 

(20141, in particular, revisited the 


previous studies assuming that the observed neutrinos could 
be produced in the inner jet of blazars and concluded that 
neither the flux nor the spectral shape suggested by the Ice- 
Cube data could be reproduced by this scenario which pre¬ 
dicts a rapid decline of the emission below 1 PeV. [Tavecchl^ 
|et al.|[^014| and [Tavecchio fc Ghisellini|[2014[ on the other 
hand, considered the distribution of lower-power blazars, 
namely, BL Lac objects and, by employing a two-zone spine- 
sheath jet model to these sources concluded that they might 
be suitable for the production of the observed PeV neutrinos 
revealed by the IceCube. 

Presently, it is very hard to define what should be the 
dominant process or the real sources that are producing the 
observed neutrino flux mainly due to the lack of more precise 
measurements. But while waiting for better measurements, 
we can explore further mechanisms and try to make reliable 
predictions in order to constrain the candidates. 

The big challenge in models that rely on hadronic pro¬ 
cesses in the AGN nuclei is how to produce the relativistic 
protons that may lead to 7 -ray emission and the accom¬ 
panying neutrino flux. Diffusive shock acceleration at the 


jet launching base was discussed by Begelman et al. 119901. 
I Levinson] ( |2000[ ) and more recently, |Vincent| ( |2014[ ) proposed 
that TeV gamma-ray emission might be produced in the BH 
magnetosphere by pulsar-like mechanisms, i.e., with par¬ 
ticles being accelerated by the electric potential difference 
settled by non uniform magnetic field. As remarked above. 


Kimura et al. (2014l discussed stochastic acceleration in an 


accreting RIAF turbulent scenario, but currently none of 
these models can be regarded as dominant or disclaimed 
given the uncertainties from the observations regarding the 
location of the gamma-ray emission (see §.5). 

In this work, we consider an alternative acceleration 
model that may occur in the vicinity of BHs which was ex¬ 
plored first in the framework of microquasars by |de Gouveia| 
Dal Pino & Lazarian 2005 (hereafter GL05) and then ex- 


tended to AGNs by |de Gouveia Dal Pino, Piovezan fc Kad-| 

owaki|26l0 (hereafter GPKIO). In this model, particles can 


be accelerated in the surrounds of the BH by the magnetic 
power extracted from events of fast magnetic reconnection 
occurring between the magnetosphere of the BH and the 
lines rising from the inner accretion disk into the corona 
(Figure 1). 

More recently, Kadowaki, de Gouveia Dal Pino & Singh 
(2015) revisited this model exploring different mechanisms 
of fast magnetic reconnection and extending the study to 
include also the gamma-ray emission of a large sample 
of sources (more than 200 sources involving blazars, non- 
blazars or LLAGNs, and galactic black hole binaries). They 


confirmed the earlier trend found by GL05 and GPKIO, ver¬ 
ifying that the fast magnetic reconnection power calculated 
as a function of the black hole (BH) mass can explain the 
observed radio and gamma-ray luminosity from nuclear out¬ 
bursts of all LLAGNs and galactic black hole binaries of their 
sam ple, spanning 10^° orde rs of magnitude in mass (see Fig. 
5 in [Kadowaki et al.|2015[ ). 

In the works above, a standard accretion disk model 


was employed, but in an accompanying work (Singh et al. 


20151, the authors adopted a magnetically dominated advec- 


tive flow (MDAF) for the accretion and obtained very similar 
results to those above, which demonstrated that the details 
of the accretion physics are not relevant in the development 
of the magnetic reconnection process and the particle accel¬ 
eration occurring in the corona. 

The magnetic reconnection acceleration model above 
(see Kadowaki et al. 20151 has been also employed to re¬ 
produce the observed multi-wavelength spectral energy dis¬ 
tribution (SED) of a few microquasars (Cygnus X-1 and 
Cygnus X-3) (Khiali, de Gouveia Dal Pino & del Valle 
(2015), henceforth KGV15) and LLAGNs (Gen A, NGG 
1275, M87 and IC310) (Khiali, de Gouveia Dal Pino & 
Sol (2015), hereafter KGS15), and the results indicate that 
hadronic mechanisms (pp and py) are the main radiative 
processes producing the observed GeV to TeV y-rays. 

The results of the works above and in particular, the re¬ 
production of the observed SEDs and the TeV gamma-ray) 
emission of the radio-galaxies by hadronic processes involv¬ 
ing particles accelerated by magnetic reconnection in the 
surrounds of the BH (KGS15), have motivated the present 
study. We aim here to calculate the spectrum of neutrinos 
arising from the interactions of accelerated protons by the 
mechanism above with the radiation and thermal-particle 
fields around the BH. According to our previous results 
(KGV15 and KGS15), these interactions produce weakly de¬ 
caying 7r° and pions. The latter may generate high energy 
neutrinos. We will then evaluate the diffuse neutrino inten¬ 
sity and compare with the IceCube data in the context of 
LLAGNs. 


For completeness, we will also compare the particle ac¬ 
celeration by magnetic reconnection with shock acceleration 
in the surrounds of the BH. 

The outline of the paper is as follows. In Section 2, we 
describe in detail our scenario. In Section 3, we describe 
the hadronic interactions and calculate the acceleration and 
radiative cooling rates. The calculation of the spectrum of 
neutrinos and their diffuse intensity for comparison with the 
IceCube data is presented in Section 4. We discuss and sum¬ 
marize our results and conclusions in Section 5. 


^ We note that the calculated reconnection power in this core 
model though large enough to explain the luminosity of galactic 
black hole binaries and LLAGNs (or non-blazars), is insufficient 
for reproducing the luminosity of most of the blazars and GRBs, 
which is compatible with the notion that the observed emission in 
these cases is produced outside the core, at the jet that points to 
the line of sight and screens any deep nuclear emission dKadowakil 
et al.|2015| |. 
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2 DESCRIPTION OF THE MODEL 

In this section we summarize the main characteristics of our 
fast magnetic reconnection model in the surrounds of the 
BH and how particles can be accelerated in the magnetic 
reconnection layer. For a more detailed description we refer 
to lKadowaki et al. 120151 


2.1 Fast magnetic reconnection in the surrounds 
of the BH 

We assume that the gamma-ray emission from low-luminous 
AGNs is produced in the core region and the particles re¬ 
sponsible for this emission are accelerated in the corona 
around the BH and accretion disk as sketched in Figure 
Turbulent dynamo inside the accretion disk or plasma 
dragging from the surrounds can build the large-scale 
poloidal magnetic fields that arise into the corona. This 
poloidal magnetic flux under the action of disk differential 
rotation gives rise to a wind that partially removes angular 
momentum from the system increasing the accretion rate 
and the ram pressure of the accreting material that will then 
press the magnetic lines in the inner disk region against the 
lines anchored into the BH horizon allowing them to recon¬ 
nect (see Fig. 1). The magnetic field intensity in the inner 
region of the accretion disk corona is approximately given 
by the balance between the magnetic pressure of the BH 


magnetosphere and the accretion ram pressure (Kadowaki 
et al.|2015| ): 

( 1 ) 


B 9.96 X ® G. 


Where rx = Rx/Rs is the inner radius of the accretion 
disk in units of the BH Schwartzchild radius (Rs) (in our 
calculations we assume rx = 6); ^ is the mass accretion 
disk rate in units of the Eddington rate (^ = M/MEdd), 
with MEdd = 1-45 X 10^®m g s“^), which we assume to be 
^ ~ m is the BH mass in units of solar mass. 

The presence of embedded turbulence in the nearly 
collisional MHD coronal flow of the core region of the 
AGNs can make reconnection very fast with a rate Vr ~ 
VAihnj/Lf^'^{vturb/vA)^, where li„j and Vturb are the injec- 
tion scale and velocity of the turbulence, respectively (e.g., 


Lazarian & Vishniac 


1999 


hereafter LV9t|_| This relation 
shows that the reconnection rate is of the order of the Alfven 
speed va , which in the systems here considered may be near 
the light speed. 

The magnetic reconnection power released by turbulent 
driven fast reconnection in the magnetic discontinuity region 
(Figureis given by (Kadowaki et al.]|2015l: 


Tl/ 1 cc in35 ,-0.5 -0.62,-0.25 7 -2,.0.75 -1 

W ~ 1.66 X 10 tp rx I Ixq 4 rn erg s , (2) 



Figure 1. Scheme of magnetic reconnection between the lines 
arising from the accretion disk and the lines anchored into the 
BH horizon. Reconnection is made fast by the presence of em¬ 
bedded turbulence in the reconnection (neutral) zone (see text 
for more details). Particle acceleration may occur in the mag¬ 
netic reconnection zone by a first-order Fermi process (adapted 
from GL05). 


where I = L/Rs is the height of the corona in units of Rs', 
lx = Lx/Rs, Lx is the extension of the magnetic reconnec¬ 
tion zone (as shown in Figure[^, g = [1 — (3/rx)°'^]°'^® and 
Va = VAoi^ is the relativistic form of the Alfven velocity, 
with vao = Rt■> R being the local magnetic field, 
p ~ ncrrip the fluid density, and ip = [1 + in 

this work, we find that vaq ~ c. 

The acceleration region in our model corresponds to 
the cylindrical shell around the BH where magnetic recon¬ 
nection takes place, as shown in Figure This shell has 
a length lx, with inner and outer radii given by Rx and 
Rx + ARx respectively, where ARx is the width of the 
current sheet given by ( Kadowaki et al.]|2015 |: 

AO ~ o 0.1 in4 ,-0.31 0.48,-0.15, -0.75,--0.15 /o\ 

ARx = 2.34 X 10 Ip rx I Ixq ? rn cm. (3) 

This magnetic reconnection power (Eq.|2| will both heat 
the surrounding gas and accelerate particles. As in |Kad-| 
|owaki et al.| ( |2015[ ), we assume that approximately 50% of 
the reconnection power goes to accelerate the particles (see 
§. 2.2.). This is consistent with plasma laboratory exper¬ 
iments of reconnection acceleration ([Yamada et al. |2014[) 


and also with the observations of solar flares (e.g., Lin fc 
Hudsun|1971 1. We further assume that this power is equally 


shared between the protons and electrons/positrons, so that 
the proton luminosity will be 25% of the calculated value by 

Eq.[2] 

The particle density in the coronal region in the sur¬ 
rounds of the BH is ( Kadowaki et al.]|2015 |: 

~or»Ov, inis -0.375/0.57-0.75 -2^0.25_-l _-3 

Tic = 8.02 X 10 ip I Q ^ ^ cm . (4) 


^ See Fig. 5 in |Kadowaki et 5I] l |201^ . Accretion rates ^ between 
0.05 < ^ < 1 are able to produce magnetic reconnection power 
values which are enough to probe the observed gamma-ray lumi¬ 
nosities from LLAGNs. 

According to the LV99 model, even weak embedded turbulence 
causes the wandering of the magnetic field lines which allows for 
many independent patches to reconnect simultaneously making 
the global reconnection rate large (for more details see [Kadowaki] 
jet al.|2015| |. This theory has been confirmed numerically by means 


of 3D MHD simulations ( Kowal et al.|2009 2012 i. 


2.2 Particle acceleration due to the magnetic 
power released by fast reconnection 

In the magnetic reconnection layer (or current sheet; see Fig¬ 
ure 1) where the two converging magnetic flux tubes move to 
each other with a velocity Vr, trapped particles may bounce 
back and forth due to head-on collisions with magnetic fluc¬ 
tuations in the current sheet. As a consequence, their energy 
after a round trip may increase by < AEjE >~ VrIc, im¬ 
plying an exponential energy increase after several round 
trips. This first-order Fermi acceleration process within re- 
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connection layers was first studied by GL05 and successfully 
tested through 3D MHD simulations with test particles in¬ 
jected in current sheets with fast reconnection driven by 


turbulence (Kowal et al. 

2011 

2012 see also ( 

de Gouveia 

Dal Pino, Kowal & Lazarian||2014||de Gouveia Dal Pino & 

Kowal||2015 

for reviews14 


From the results of the 3D MHD numerical simulations 
(Kowal et al. 2012), we find that the acceleration rate for a 
proton is given by (see also KGV15): 


^acc,M.R.,p — 1-3 ^ lo 1 (5) 


where E is the energy of the accelerated proton, Eq = nipP, 
rrip is the proton rest mass, to = laccjvA is the Alfven time, 
and lacc is the length scale of the acceleration region and for 
electrons this rate is (KGV15): 


t 


-1 

acc,A'I. R. ,e 



-0.4 


t 


-1 
0 1 


( 6 ) 


where rrie is the electron rest mass. 

As stressed in GL05 (see also KGV15 and KGS15), it 
is also possible that a diffusive shock may develop in the 
surrounds of the magnetic reconnection zone due to coronal 
mass ejections released by fast reconnection along the mag¬ 
netic field lines, as observed in solar flares. A similar picture 


has been also suggested by e.g., Romero et al. (2010b i. In 


this case, the shock velocity will be predominantly parallel 
to the magnetic field lines and the acceleration rate for a 
particle of energy i? in a magnetic field B, will be approxi¬ 
mately given by (e.g., Spruit|1988 1 : 


TjecB 
E ’ 


( 7 ) 


where 0 < ?7 ^ 1 characterizes the efficiency of the accel¬ 
eration. We fix = 10“^, which is appropriate for shocks 
with velocity Vs ~ 0 . 1 c, which are commonly assumed in 
the Bohm regime (Romero et al.|2010bK 


Table 1. Three sets of model parameters for LLAGNs. 


Parameters 

Model 1 

Model 2 

Model 3 

m BH mass (Mq) 

10' 

10“ 

i(F” 

p Injection spectral index 

1.9 

1.7 

2.2 


3.1 pp collisions 


The charged pions can be created through inelastic collisions 
of the relativistic protons with nuclei of the corona that 
surrounds the BH and the accretion disk by means of the 
following reactions ( Atoyan fc Dermer|2003 Becker|[2008 1 

p + p ^ nx{T\^ + n~) + n2'K'^ + p + p ( 8 ) 

where ni and n 2 are multiplicities, tt® —>■ 7 -I- 7 (Stecker 


1970 19711, carrying 33% of the accelerated proton’s energy. 


The charged pions tt* then decay and produce neutrinos via 
7 r+ —>• and-TT” —>• Vp-\-Up-\-v^-\-e~ Vp, 

Vp, and are the muon neutrino, muon antineutrino, and 
electron neutrino, respectively ([Margolis et al.|1978[|Stecker| 


1979 Michalak et al. 19901. The pp cooling rate is almost 


,-i 

^pp — TLiCCJp 


independent of the proton energy and is given by (Kelner 

[ 200 ^ 

( 9 ) 

where, rii is the coronal number density of protons which 
can be calculated by eq. and kpp is the total inelasticity of 
the process of value ~ 0.5. The corresponding cross section 
for inelastic pp interactions Upp can be approximately by 
( Kelner et al.|2009 l 

2 


app{Ep) = (34.3 -b 1.88Q + 0.25Q^ 


1 - 


Ep 


mb, ( 10 ) 


where mb stands for milli-barn, Q = In ^ iTev ) ’ 

ton threshold kinetic energy for neutral pion (P) production 
is Eth = 2mpP(l + j^) ~ 280 MeV, where rriTTp = 134.97 
MeV is the rest energy of 7 r° ( Villa fc Aharonian|2009 1. This 
particle decays in two photons with a probability of 98.8%. 


3 HADRONIC INTERACTIONS 

In KGS15, we have demonstrated that the core region of 
LLAGNs is able to accelerate protons up to energies of a 
few lO^^eV through the first-order Fermi magnetic recon¬ 
nection mechanism described in the previous section. This 
indicates that these sources could be powerful CR acceler¬ 
ators. We show below that these protons can cool very effi¬ 
ciently via synchrotron, p'y and pp interactions in the region 
that surrounds the BH of these sources (Figure]^. As re¬ 
marked, these hadronic interactions lead to the production 
of HE 7 -rays and HE neutrinos via decays of neutral and 
charged pions, respectively. In KGS15, we have calculated 
the spectral energy distribution of the HE 7 -ray emission 
for the LLAGNs for which this emission has been detected. 
Below, we calculate the HE neutrino emission from the nu¬ 
clear region of these sources. 

® Particle acceleration within reconnection sheets has been also 
extensively tested numerically in collisionless fluids by means 
of 2D (e.g., [Zenitani &: Hoshin^|2001[ [Drake et al.||2006| |2010[ 
Cerutti et al.|2013| l and 3D particle in cell simulations ( jSironi fej 

Spitkovsky|2014 1 . 


3.2 p 7 interactions 


The photomeson (py) production takes place for photon en¬ 
ergies greater than Eth ~ 145MeV. Pions are also obtained 
from the P 7 interaction near the threshold via the channels 
I Atoyan fc Dermer|2003' I 

p-b 7 —t p-b 7r°, (11) 

with 7 r° —>■ 7 -b 7 carrying 20 % of accelerated protons energy 
and 


p-b7—tp-b-TT^-bTr' 


( 12 ) 


and the charged pions will also decay producing neutrinos 
as described in §. 3.1. 

The radiative cooling rate for this mechanism in an 
isotropic photon field with density nph{Eph) can be calcu¬ 
lated by ( |Steckerj|1968| ): 


tpj (Ep) 


27| 


L 


M 

th 

STp 


dE, 


'ph - 


'^ph 


■ph) 


Elh 


X / derO-^7 (erjK^^ (erjcr, 


(13) 
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Figure 2. Acceleration and cooling rates for protons in the core regions of LLAGNs with a central black hole mass (a) M = IO^Mq 
(M odel 1), (b) M = 10®M© (Model 2), and (c) M = 10®M© (Model 3). 


where in our model the appropriate photons come from the 
synchrotron radiation n '^phi^Eph'j — n 

synchi^^') ^ ^^.2 1 

er is the photon energy in the rest frame of the proton, and 


is the inelasticity of the interaction. Atoyan & Dermer 


( |2003[ ) proposed a simplified approach to calculate the cross- 
section and the inelasticity which are given by 


, N ^ / 340 ^barn SOOMeV < < 500MeV 

cTp^ier) « I ^20 ^barn tr > SOOMeV, 

and 

, , / 0.2 SOOMeV < Cr < SOOMeV 

Kp^{er) « j Q g ^ SOOMeV. 


(14) 


(15) 


® We find that for photomeson production, the radiation from 
the accretion disk is irrelevant compared to the contribution from 
the coronal synchrotron emission above the disk (see KGV15 and 
KGS15 for a detailed derivation of the synchrotron rate and its 
radiation field density). 


4 NEUTRINO EMISSION AND DIFFUSE 
INTENSITY 


To calculate the neutrino emission from the nuclear region 
of an LLAGN, we consider a population of protons acceler¬ 
ated by magnetic reconnection in the surrounds of the BH 
according to the model described in Section 2. 

We assume for these accelerated particles an isotropic 
power law spectrum (in units of erg“^cm“®s“^) (see e.g. 
KGV15 and references there in): 


Q{E) = QoE ^exp[-E/Emax] 


(16) 


where p > 0 and Emax is the cut-off energy. 

The normalization constant Qo above is calculated from 
the total power injected to accelerate the protons according 
to the relation: 


l^rf- 

Jv j 


dE E Q{E) 


(17) 


where V is the volume of the emission region around the 
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magnetic reconnection zone and Lp corresponds to the mag¬ 
netic reconnection power W given by Eq. To calculate W 
we have adopted the following snitable set of parameters 
^ = 0.7, Rx = 6Rs, Lx = 10J?s, and L = 20i?s. 

The maximum energy of the accelerated particles Emax 
is derived from the balance between the magnetic reconnec¬ 
tion acceleration rate (Eq. 1^ and the radiative loss rates as 
given in Section 3. Figure [Compares these rates for protons 
considering LLAGNs with three different BH masses 10^, 
10® and 10® Mq. We have also considered different power- 
law indices (p) for the injected particle spectrum in each of 
these models (see Table 1) which are compatible with the 
values derived from analytical and numerical studies of first- 
order Fermi acceleration by magnetic reconnection and also 
with values inferred from the observations (e.g., KGS15). 
The calculated values of 73, W, ARx and ric from Eqs. 
l^for these three representative source models are listed in 
Table For simplicity, we consider the derived proton lu¬ 
minosities (which are ~ 1/4W) and the emission properties 
of these three models to characterize the whole range of 
LLAGNs in the calculation of the HE neutrino flux below. 
The adoption of this approach, rather than accounting for 
a whole range of BH mass sources allows us to avoid the 
introduction of further free parameters in the modelling. 

In Fig. for comparison we have also calculated the 
proton acceleration rate due to a shock formed in the sur¬ 
rounds of the reconnection region (Eq. for the same set of 
parameters as above. As in KGV15 and KGS15, we find that 
the maximum energy attained from magnetic reconnection 
acceleration is higher than that from the shock . It should be 
also remarked that protons with these calculated maximum 
energies have Larmor radii smaller than the thickness of the 
reconnection layer ARx (eq. [^, as required. 

The neutrinos that are produced from pion decay will 
escape from the emission region without any absorption and 
their spectrum is given by (Tavecchio et al.|2014 Kimura et 

aL][MI| ): 


E,L^{E,) ~ (0.5t-p' -b 


EpI/p , 


(18) 


where E^, is the neutrino energy and Ep the proton energy. 
Since Figure [^demonstrates that the pq emission cools the 
protons faster than pp collisions, the dominant hadronic pro¬ 
cess in our model is the py emission. Therefore, this mecha¬ 
nism will prevail in the production of the neutrinos and the 
first term of eq.j^can be neglected. In py interactions, Ei, is 
related with the parent proton energy through the equation 
Eu — 0.05Ep (Spurio 20151, because the average energy of 


the pion is ~ 0.2 of the parent proton energy and in the de¬ 
cay of the tt'*' chain four leptons are produced (including one 
electron neutrino as remarked), each of which has roughly 
1/4 of the pion energy. It has been also demonstrated in 
[s5 urio| ( |2015[ ) that the ratio of the neutrino luminosity to 
the photon luminosity from py interactions is ~ 1/3. 

In consistency with the statement above, the maximum 
energy of the produced neutrinos can be calculated from 
Eu.max = 0.05Er,.max l|Becker||2008l iHazlen 20071, which 


according to our model is ~ 3 x lO^^’eV for a source with a 
black hole mass Mbh = lO^M©, ~ 5 x 10^®eV for a source 
with Mbh = 10® Mq, and ~ 2 x 10^®eV for a source with 
Mbh ~ 10® Mq. 

The total diffuse neutrino intensity from the extragalac- 



Figure 3. Calculated diffuse intensity of neutrinos from the cores 
of LLAGNs considering our magnetic reconnection acceleration 
model to produce the protons and gamma-ray photons for three 
different BH masses. The data are taken from IceCube measure¬ 
ments ([Aartsen et al.|2014^. 


tic sources we are considering here, i.e., LLAGNs may have 
contributions from different redshifts. Neglecting evolution¬ 
ary effects in the core region of these sources, we can esti¬ 
mate the total intensity as iMurase et al.|2014| 


-!>. = 


dz 


f 

J L. 


'\/(l + z)^Q,m + 


AtvHq Jq 

Lmax 

X / dL-y py \E~f 5 Zj ~ 

Eu 


(19) 


where L-y is the y-ray luminosity, and pj{Ly, z) is the y-ray 
luminosity function (GLF) of the core of the sources, de¬ 
fined as the number density of sources per unit comoving 
volume, per unit logarithmic luminosity between the red- 
shifts 2 = 0 to 2 = Zmax, being the latter the maximum 
observed redshift for radiogalaxies, Zmax 
|al.||2005]). GLF is integrated from L. 


~ 5.2 (Klamer et 
^rmn to Ljyiax which are 
obtained from Fermi-hAT observations and are given by 


10 and 10 erg/s, respectively (Di Mauro et al. 20141. 
The values for the cosmological parameters are assumed as: 
Ho = 70 km s“'^Mpc“'^, = 0.3 and Ha = 0.7. 

We evaluate the GLF, py(Ly,z) as in Di Mauro et al.j 
( |2014[ ), from the estimated radio luminosity function (RLF) 
which for non-blazars is given by 

P^{Ly,z) = Pr,tot{Lr':tot\Ll%^ULy)),z) 


d log I/i 


5GHz 

r.core 


dlogL®?„t 


( 20 ) 

dlogL^ dlogZ/SGHij,' 

dlogL®™/*e/fHogI /7 and dlog Z/®|/,S^/dlog can be 

calculated by ( Di Mauro et al.|2014[ ) 

logL^ = 2.00 ± 0.98 -b (1.008 ± 0.025) log (L®™//), (21) 

and 

logL®™// = 4.2 ± 2.1 -b (0.77 ± 0.08) log (L®;?,?), (22) 

where L/ 
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Table 2. Physical conditions around the LLAGNs represented by models 1, 2 and 3, 
obtained from Eqs. 1 to 4, using = d, I = 20, lx = 10 and g = 0.7. 



Parameters 

Model 1 

Model 2 

Model 3 

B 

Magnetic field (G) 

2.8 X 10'‘ 

8874 

2806 

W 

Magnetic reconnection power (erg/s) 

2.4 X 10"*® 

2.4 X 10"*® 

2.4 X 10"*"* 

ARx 

Width of the current sheet (cm) 

7.2 X IQi® 

7.2 X 10®® 

7.2 X 10®"* 

TLc 

Coronal particle number density (cm~^) 

3.6 X IQi® 

3.6 X 10® 

3.6 X 10® 


core luminosities, respectively. The total RLF, 
z), is found from interpola¬ 
tion of the observed data for radio-galaxies provided by 
[Yuan &: Wang| ( f2012[ ): 

PrMLr'forMLj), z) = (-1.1526 ± 0.0411) logL^™', 

-h (0.5947 ± 0.1224)2 -h 23.2943 
± 1.0558 Mpc-"(logL®™l)”'. 

(23) 

The resulting neutrino flux is shown in Figure It 
was calculated using eq. above, considering the maxi¬ 
mum neutrino energies obtained for sources with the three 
different BH masses (as in Figure 2). 

Sources with Mbh = iGM q result a spectrum that 
matches better with the observed most energetic part of the 
neutrino flux by the IceCube, at ~ 3 x 10^® to 10^® eV, while 
sources with BH masses of the order of 10*Mq produce a 
spectrum that nearly fits the observed neutrinos flux in the 
range of lO^'^ — 10^®eV, and sources with mass ~ 10 ®Mq 
the narrow energy band 5 x lO^^eV — 10^‘^eV as well as the 
upper limit at 5 x lO^^eV. 


5 DISCUSSION AND CONCLUSIONS 


In this work we have explored a model to describe the ob¬ 
served flux of extragalactic very high energy (VHE) neutri¬ 
nos by the IceCube ( Aartsen et al.|[2014 l in the framework 
of low luminosity AGNs (LLAGNs), or more specifically, of 
radio-galaxies. The recent detection of gamma-ray emission 
in the TeV range in these sources makes them also poten¬ 
tial candidates of VHE neutrino emission via the decay of 
charged pions which can be produced by the interaction of 
accelerated relativistic protons with ambient lower energy 
photons and protons. 

We have examined here a fast magnetic reconnection 
mechanism in the surrounds of the central BH occurring be¬ 
tween the lines lifting from the accretion disk into the corona 
and those of the BH magnetosphere to accelerate particles 
to relativistic energies through a first-order Fermi process in 
the reconnection layer (GL05, Kowal et al.|2012 l. Recently, 
it has been demonstrated that this model successfully re¬ 
produces the observed gamma-ray luminosity of hundreds of 


LLAGNs (Kadowaki et al. 2015 and SGK15) and also shapes 


the SEDs of several radio-galaxies, particularly reproducing 
their TeV gamma-ray energies mainly via photomeson (py) 
interactions (KGS15). 

Applying the same acceleration model as above (see 
Section 2), considering three different BH masses, we have 
shown that also the observed VHE neutrino Icecube flux 


(Aartsen et al. 20141 can be obtained from the decay of 


charged pions produced in photomeson interactions involv¬ 
ing the accelerated protons and Synchrotron photons in the 
core region of these sources (Figs. 2 and 3). 

Specifically, in Fig. we compared the magnetic re¬ 
connection acceleration rate (derived from the numerical 
simulations of [Kowal et al.| ( [2012 1 and calculated for the 
source parameters) with the relevant hadronic cooling pro¬ 
cesses and obtained the maximum energy for the accelerated 
protons mainly constrained by the py interactions. In Fig.[^ 
we also compared the magnetic reconnection with the shock 
acceleration rate in the surrounds of the BH for the same 
parametric space and demonstrated the higher efficiency of 
the first process in this region. According to our results in 
Fig.m protons are able to accelerate up to energies of the 
order of ~ lO^^eV and therefore, are suitable to produce 
0.1-1 PeV neutrinos. 

Fig. [3] indicates that the observed neutrino flux in the 
few PeV range can be matched by sources with Mbh ~ 
iGM q (Model 1), while the flux in the energy range of 
0.1PeV< Eu <lPeV can be matched by sources with 
Mbh ~ 10® Mq (model 2), and that in the range < O.lPeV 
can be fitted by sources with Mbh ~ 10®eV (model 3). 

We note that, although the calculated neutrino flux 
was obtained from the integration of the contributions of 
LLAGNs over the redshifts between z=0 and 5.2 (eq. |19| ) 
considering, for simplicity, sources with only three charac¬ 
teristic values of BH masses, one may naturally expect that a 
continuous integration considering the sources with all possi¬ 
ble BH masses within the range 10^ —10®M© should provide 
a similar fitting to the observed data. We also note that our 
model is unable to explain the IceCube upper limits at the 
~ lOPeV range (also depicted in Fig. [^, which are probably 
due to other astrophysical compact source population. 

Furthermore, we expect that with the 10-fold increased 
sensitivity at TeV energies, and the larger field of view and 
improved angular resolution of the forthcoming gamma-ray 
observatory CTA ( jActis et al. 112011] [Acharya et al.|[2013p , 
the list of LLAGNs with confirmed detection of gamma- 
ray emission at TeV energies (which currently has only four 
sources: Gen A, Per A, M87 an IC310), will increase sub¬ 
stantially, allowing for a more precise evaluation of the con¬ 
tribution of individual sources for the IceCube neutrino flux. 

As remarked in §. 1, other models have been proposed 
in the literature to explain the IceCube neutrino flux which 
cannot be discarded or confirmed, considering the current 
poorness of the data available. 


Tavecchio et al. (20141 and Tavecchio & Ghisellini 


(20141, for instance, have proposed that the lower power 
blazar class of BL Lac objects could be promising candi¬ 
dates to produce the observed neutrino flux. In their two- 
zone jet model, the neutrinos are produced by photomeson 
interactions involving photons emitted in the slower, outer 
layer that envelopes the faster inner jet component. A lim- 
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itation of this model is that the high-energy cut-off of the 
accelerated protons, as well as their injected power are free 
parameters, unlike in our model where both quantities are 
directly obtained from the magnetic reconnection accelera¬ 
tion mechanism. Besides, since the BL Lacs are a subclass 
of the blazars, another difficulty with this model is that it 
is not clear whether the remaining more powerful blazars, 
which are also TeV gamma-ray emitters, can or cannot pro¬ 
duce neutrinos. According to the recent studies of |Dermer et| 
al. (2014 1 and Murase et al. (20141, which employed a single 


zone jet model, the powerful blazars would not be suitable 
candidates to explain the IceCube data. These analyses and 
the relatively large number of free parameters employed in 
the evaluation of the neutrino flux leave the question on 
whether or not blazars do contribute to the IceCube data 
opened. 


Another model to explain the observed neutrino flux 


has been proposed by Kalashev et al. (20141 who studied 


photo-pion production on the anisotropic photon field of a 
Shakura-Sunyaev accretion disk in the vicinity of the BH 
in AGNs. But this model does not provide an acceleration 
mechanism either and therefore, the proton high energy cut¬ 
off is also a free parameter. 


Recently, radio galaxies have been also discussed as pos¬ 
sible sources of the observed HE neutrinos by Becker et al. 
(2014). They demonstrated that FR I radio galaxies would 
be more probable sources of this emission than FR II radio 
galaxies. In this work, as we considered the global diffuse 
contribution from LLAGNs spread over a range of 2 values, 
we cannot distinguish the relative contributions from both 
classes. 


Finally, another recent study ( Kimura et al. [2(514 1 spec¬ 
ulates that the protons responsible for the neutrino emis¬ 
sion could be accelerated stochastically by the turbulence 
induced in a RIAF accretion disk in the core region of 
LLAGNs. This acceleration process should be essentially a 
second-order Fermi process and therefore, less efficient than 
a first-order Fermi process. Nevertheless, their analytically 
estimated acceleration rate face « £(“ 0.35 gggjjjg i^g 
large when compared to that predicted for first-order Fermi 
processes, as for instance in the present study where the ac¬ 
celeration rate has been extracted directly fr om 3D MHD 
simulations with test particles oc E 


Kowal et al. 


M 2 I KGV15), or in shock accelerati on (fo r which analytic 


predictions give tjj-c oc E~^ (Spruit 1988 1 ). Furthermore, 
since the candidates to produce PeV neutrinos in this case 
are radio-galaxies, which are the observed 7 -ray emitters, 
it seems that the employment of the gamma-ray luminos¬ 
ity function (GLF; [Di Mauro et al. 20141 to calculate the 
diffuse neutrino intensity as we did here seems to be more 
appropriate than the employment of the luminosity function 
in X-rays, as these authors considered. 


In summary, in spite of its simplicity, the numerically 
tested acceleration model applied to the core region of 
LLAGNs here presented indicates that LLAGNs are very 
promising candidates to explain the IceGube VHE neutri¬ 
nos. 
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